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ABSTRACT 
We  Investigate  the  excitation  and  propagation  of  the  three  dimensional 
electromagnetic  field  over  an  infinite  corrugated  plane  which  is  approximated 
"by  an  anisotropic  impedance  "boundary  condition.   Emphasis  is  placed  upon 
effects  of  surface  anisotropy  which  are  not  evident  in  two  dimensional  treat- 
ments.  In  particular  we  consider  the  excitation  by  a  magnetic  point  dipole 
in  detail.   It  turns  out  that  the  fields  are  determined  "by  a  scalar  wave 
function  which  satisfies  a  mixed  boundary  condition  involving  a  linear  com- 
bination of  the  wave  function,  its  normal  derivative  and  its  second  order 
tangential  derivative.   The  exact  formal  solution  is  first  derived,  and  then 
the  radiated  far  field  and  the  surface  wave  far  field  are  evaluated 
separately.   Both  the  phase  and  the  amplitude  of  the  excited  surface  wave 
are  dependent  upon  the  direction  of  obserAration.   N-umerical  results  are 
given.  The  physical  significance  of  this  solution  is  discussed.  A  com- 
parison is  made  between  this  problem  and  the  theory  of  ship  waves. 
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1.   Introduction 

In  this  paper  we  study  the  excitation  and  propagation  of  the  three 
dimensional  electromagnetic  field  over  an  infinite  corrugated  plane  which 
may  he  approximated  by  an  anisotropic  impedance  boundary  condition.   This 
investigation  is  motivated  by  its  potential  application  to  scanning  surface 
wave  antennas  .   The  fields  produced  by  a  magnetic  point  dipole  will  be 
considered  in  detail.   The  boundary  condition  is  such  that  surface  waves 
are  generated  and  it  may  be  described  mathematically  as: 

E   =ZH   =(R-iX)H    n 
X         z  z    I 

r  y  =  0 

E   =  0  J 

z 

where  the  time  dependence  is  assumed  to  be  of  the  form  e    ,   and  Z,    R,  and 

X  are  the  impedance^  resistance^  and  reactance  of  the  surface,  respectively. 

Insofar  as  siirface  waves  are  concerned,  the  most  important  case  is  such  that 

R  and  X  are  positive  in  sign  and  R  is  much  smaller  than  X,  i.e.,  Z  is  almost 

purely  reactive  . 

It  should  be  noticed  that  any  impedance  boiuidary  condition  approximating 

a  dielectric  coated  surface  is  isotropic.   The  electromagnetic  fields  excited 

by  a  uniform  line  source  above  an  infinite  reactive  plane  have  been  investigated 

2  3 

by  Cullen  ,  and  Friedman  and  Williams  .   The  propagation  of  oblique  surface 

waves  over  a  corrugated  plane  has  been  considered  in  the  literature,  but  it 
has  not  been  treated  three-dimensionally  .   The  solution  of  our  three  di- 
mensional problem  will  show  the  effect  of  the  surface  anisotropy  which  is  not 
evident  in  the  two  dimensional  problems.   The  problem  we  treat  may  be  formulated 
by  the  use  of  a  double  Fourier  transform.   However,  we  reduce  the  problem  to 


one  involving  only  a  single  Fourier  transform  "by  introducing  an  auxiliary 

function  which  is  similar  to  those  recently  employed  for  several  other 

4567 
problems  in  diffraction  theory  '    '    '    . 

In  section  2  we  obtain  the  formal  solution  for  the  fields  of  a 
magnetic  point  dipole  oriented  in  the  direction  of  corrugations  and  located 
above  a  coringated  plane  which  is  approximated  by  an  anisotropic  impedance 
boundary  condition.   The  dipole  is  first  decomposed  into  phased  line  sources, 
and  the  fields  of  the  line  sources  are  then  synthesized  into  the  dipole 
fields.   In  section  3  the  far-zone  radiation  field  is  evaluated  by  an  algebraic 
method  involving  asymptotic  differentiation  and  the  excited  surface  wave  far  field 
is  determined  by  applying  the  method  of  stationary  phase.  All  the  surface 
wave  characteristics  which  include  the  amplitude,  the  phase,  and  the  ex- 
ponential decay  factor  are  functions  of  the  direction  of  observation  and  the 
surface  reactance.   In  section  k   we  summarize  the  peculiar  properties  of  the 
three  dimensional  surface  wave,  and  present  the  svirface  wave  power  patterns 
and  the  surface  wave  phase  patterns.   A  comparison  is  also  made  between  this 
problem  and  the  theory  of  ship  waves. 
2.   The  Formal  Solution 

In  this  section  we  proceed  to  study  the  electromagnetic  field  which 
arises  from  a  magnetic  point  dipole  lying  above  the  surface  of  a  plane 
characterized  by  an  anisotropic  impedance  boundary  condition.   Consider  a 
magnetic  point  dipole  located  at  the  point  X=0,  y=h,  z=0  (See  Fig.  l) 
and  oriented  in  the  direction  of  the  z-axis.   We  are  interested  in  the  free 
space  region  y  >  0.   The  impedance  boundary  condition  prescribed  on  the  plane 
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siirface  y  =  0  is  given  by 


(2.1) 


E        =     0 
z 


E     =     ZH     =-iXH 
X  z  z 


where  the  impedance  Z  is  supposed  to  be  piirely  inductive  and  will  support 
surface  waves.  We  now  wish  to  solve  the  time  reduced  Maxwell's  equations 
subject  to  the  prescribed  singularity  and  boundary  conditions. 

If  we  assiMie  the  time  dependence  to  be  of  the  form  e    ,   the  mono- 
chromatic Maxwell's  equations  maybe  reduced  to: 


(2.2) 


curl  H 


—  iwe  E 


curl  E  =  i  to  M  H  -  5(x)  5(y  -  h)  S(z)  z 


where  e  and  f-i  are  the  permittivity  and  magnetic  permeability  of  free  space. 
Because  of  the  boundary  condition,  the  z-component  of  the  electric  vector 
is  absent  in  the  fields  and  hence  the  field  is  completely  determined  by  the 


value  of  a  scalar  wave  function  f .   We  have 

1 


(2.3) 


E  =  curl  f  z  and  H  =  - 


iu)  n 


curl  curl  f  z 


I.e. 


(2.U) 


and 


(2.5) 


E 


H 


H 


H 


0,   E  =  ^   ,        E 

^      ay    ^      y 


a^f 


ioj(i  SxSz 

1      a2^ 


iU)  [i  SySz 


1 

iwH 


,-az2 


af 

dx 


+  k  f 


The  function  f  satisfies  the  equation 


.1^- 


(2.6)  (V^  +k^)f  =  5(x)   5(y  -  h)  5(z) 

where  v  is  the  three-dimensional  rectangular  Laplacian,  k  is  the 
propagation  constant  of  free  space,  and  S  is  the  Dirac  Delta  function. 

The  appropriate  boundary  condition  maybe  deduced  from  (2.1)  and  (2.3): 

2 

(2.7)  1^  +  -  (^2  f  +|4)  =  0   when  y  =  0 

where  —  is  a  positive  number.   The  mathematical  problem  has  been  reduced 
to  that  of  solving  the  inhomogeneous  wave  equation  (2.6)  subject  to  the 
mixed  second  order  boundary  condition  given  by  (2.7).   Furthermore,  the  far 
field  should  satisfy  the  radiation  condition. 

Applying  Fourier  transform  to  eqs.  (2.6)  and  (2.7),  we  obtain 

2      2 

(2.8)  (-^+  -^  +k2)F  =  S(x)  &(y-h) 

Sx     ^ 

(2.9)  ^  +  ^  K^  F  =  0 

where  the  functions  f(x,y, z)  and  F(x,y,s)  are  related  by  the  following  Fourier 
transform  pair: 

00 

(2.10)  f(x,y,z)  =  —   /    F(x,y,s)   e^  ^  ^  ds 

—  00 
00 

(2.11)  F(x,y,s)  =  j        f(x,y,z)  e""-  ^  ^  d  z 

-00 

2    2    2 
and  K  =  k  -  s  .   Mathematically  the  transformed  equations  are  a  two- 
dimensional  inhomogeneous  wave  equation  and  a  mixed  first-order  boundary  con- 
dition.  Physically  the  point  dipole  has  been  decomposed  into  phased  line 


•5- 


sources.   We  are  going  to  first  solve  the  problem  of  a  phased  line  soiorce 
and  then  synthesize  the  line  source  fields  into  the  dipole  field. 

The  field  of  a  phased  line  source  may  be  analyzed  either  by  the 
Fourier  transform  method  or  by  using  another  auxiliary  function  which  is  a 
linear  combination  of  the  unknown  wave  function  and  its  derivative.   We 
choose  to  employ  the   latter  method  because  it  gives  another  example  of  this 
relatively  new  technique.   Let  us  introduce  the  auxiliary  function^ 

(2.12)  G  =   (—  +  X)   F 

X   2 
where  \   =     —  K  .   Then  G  satisfies  the  equation 

2       2 

(2.13)  (^  +  -^+k2)G  =  -  (-^  -  X)  &(x)  &(y-h) 

bx.  by  bh 

and  the  simple  homogeneous  boundary  condition 

(2.14)  G  =  0  when  y  =  0 

2     2 

We  first  treat  the  case  of  positive  X^i.e.^s  <  k.   If  the  function 

F  is  denoted  by  F  when  X  is  positive,  then  F-,  can  be  found  from  G  by  integrating 
(2.12)  in  the  following  manner. 

y 

(2.15)  \U,Y,x)     =     e"^  y  jT  e^^  G(x,ri,s)  dy  +e-^y  T(x,s) 

—  00 

where  T(x,s)  is  an  arbitrary  function  of  x  and  S,  which  except  for  certain 
discontinuities  in  the  derivative  must  be  a  homogeneous  solution  of  (2.13 ). 
The  discontinuities  must  be  adjusted  so  that  F^(x,y,s)  meets  the  requirements 
of  the  problem.   This  adjustment  follows  after  eq.  (2.l8)  below. 
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S 


ion 


By  applying  the  operator  -  (^^ "  ^)  "to  the  unique  Green's  funct 
satisfying  eq.  (2,8)  and  vanishing  at  y  =  0,  we  obtain  the  solution  for 
the  auxiliary  function: 


(2.16)  a  =  i(^-x){  H 


(1) 


K\/x2+  (y-h)2 


-H 


(1) 


K\/x2  +  (yfh)2 


We  substitute  (2,l6)  into  (2.15)  and  introduce  the  excited  surface 
wave  terms  to  play  the  role  of  e    T(x;s) 


y 


(2.17)  v--^)  = -i^"!/ ^'N^^o 


(i)r 


V^s^H-h)^ 


+  ^  H 

an   o 


K/x2+(r+h  )2 


dTl 


y 


7 


^^Xe-^M'    e^^ 


(1) 


"  ^ I    (1)  r  ^ 


(Tl+h)2J     M^ 


+     C  e-^y  +  iVK2  +^2 


X 


Integrating  by  parts,  we  have 


i    xx(l) 


(2.18)      F^(x,y,s)   =  -  I  H^^^     ^x2   +  (y_h)2 


i|       o 


K 


yx2    +  (y  +  h)^ 


2  J  o 


KVx2   +  (n   +h)2 


d  Tl 


-,     c     e-^^"^  Vk^   +x2      |x| 
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Now  the  siirface  wave  amplitude  C  may  be  determined  by  requiring  the 

continuity  of  — —  across  the  line  x  =  0^  i.e.^ 
ox 


(2.19) 


=   Lim  \      —  F   (x  >  0,  y)  -  ^  F  (x  <  0,  y) 
x^O  I     dx    -^  dx  0. 


=  0 


Substituting   (2.18)   into  the   jump  condition   (2.I9)  yields 


(2.20) 


Lim 
x-»0 


^^K^"^''/    ^^M""^    [Vx2   +(Tl+h)2]    x/x2   +  (11   +  h)2 


dTl 


+  i   2  C^j^   _^^2     e"^y    ■"  ^\/k2    +X^    \x\ 


■      =      0 


The  limit   of  the  integral  in  the  above  equation  contributes  only  in  the 
neighborhood  of  T]     =  -  h.      Hence 


(2.21) 


X  K  e 


-Xh 


-h  +e 


Sy  k2   +  X 


Lim 
x^O 


H. 


(1) 


[Kyx^  +  (n  +hf 


X 


-h-e 


yx^   +  (Tl   +  h)2 


dTl 


Using  the  asymptotic   fomiula  of  Hahkel  function  for  small  argument,  we   find 
that 


(2.22)       C 


X  K  e 


•h   +  e 
■Xh                    p  _2i  xdTl 

Lim     /        ( -) 


r  -   'd  X 

)J         JtK/x2   +  (T1+] 


2yK2   +  x2       x-*oJ         rtl^x2   +  (Tl+h)2  Vx^   +  (11+  h)^ 

-h-e 


(2.23) 


Carrying  out  the  integration  gives 

.  ^     -Xh 
1  X  e 


VK^   +  x2 


2    2 
Nov  we  consider  the  case  of  negative  X,  i.e.,  s  >  k  .   If  the 

function  F  is  denoted  by  F  when  X  is  negative,  then  F  may  be  found  from 

G  by  integrating  (2.12)  as  follows: 

oo 

(2.24)  '   F2(x,y,s)  =  -  e'^^  f      e^^  G  (x,Tl,s)  dTl 

y 

Substituting   (2.l6)   into   (2.2U)  and  integrating  by  parts,  we  have 


F2(x,y,s) 

1   (1' 

.   (1) 

" 

(2.25) 

00 

) 

4      o 

L^n/x^    +   (y 

.h)2j 

- 1  ^  -'"■'  J  -'■"  \ 

_K    s/^ 

+  (in  +h)2 

dTl 

y 

The  limits  of  integration  in  (2.24)  have  been  so  chosen  that  they 

r^F 
insure  the  convergence  of  the  integral.   Because  §_  is  continuous  across  the 

ax 

the  line  x  =  0,  no  siorface  waves  are  excited  in  this  case.   Since  the  arguments 
of  Hankel  Functions  in  (2.28)  are  imaginary,  F  will  not  give  rise  to  a 
radiation  field. 

The  substitution  of  (2.18)  and  (2.25)  into  (2.10)  completes  the  formal 
solution  of  this  problem.  It  is  difficxxlt  to  evaluate  the  integrals  exactly; 
however,  we  shall  obtain  the  far  field  expressions  in  the  next  section. 

3.   Determination  of  the  Far   Fields 

In  this  section  we  determine  the  radiated  far  field  and  the  surface 
wave  far  field.   We  first  treat  the  radiated  far  field  by  an  algebraic  method 
involving  asymptotic  differentiation.   This  method  has  proved  useful  in  two 
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dimensional  problems  '    '        .      It  gives  the  radiated  far  field  directly, 


without  intermediary  use  of  the  previously  obtained  complete  solution. 
Let  us  now  introduce  the  following  function 


(3.1) 


dy    w  M.        az2 


Then  g  satisfies  the  eq_uation 
(3.2)        (v^   +k2)  g     = 


S  X         o         S^     - 

—     +  —   (k^   +  ^— ^) 
dh         u)  U  dz  ^  . 


5(x)   &(y-h)   6(z) 


and  vanishes  when  y  =  0.      The   solution  for  g  may  be  written  down  immediately: 


(3.3)        g    =    r" 


>.         r        n         :^2    -[    p    ik  R' 
Sh      oJM.  Sz      -I     L      r' 


ikR'S 


where   R'   =  \/x^  +  (y  -  h)^+  z^  fn     R  -  h  sin  0  sin  (j)  , 

R'  •  =\/  x2  +  (y  +  h)^  +  z^  ~  R  +  h  sin  9   sin  (j)  , 
R   =  \/x2  +  y2  +  z^   ,  9   =   tan" 


^         -1   y 

,   and  (p  =  tan         — 
x'"+  y""  ^ 


y3775 


We  may  expect  the  far-zone  radiation  field  of  f  to  be  of  the  form 
(3.^) 


mje^)   ikR 
R 


Substituting  (3.3)  and  (3.^)  into  (3.1), introducing  polar  coordinates,  and 
keeping  terms  of  order  -x,   we  obtain 

fr  -1  ■      a      ■      ^       X  /  ?  ?    2  flj  -ikh  sin  0  sin  (t) 
J  -ik  sm  9   sm  <\)  -   —  (k  -k  cos  0)  e  ■ 


(3.5) 


X    2   2    2 
-ik  sin  0  sin  d)  +  -— ,  (k  -k  cos  0) 

'     u)  U 


^ikh  sin  9   sin  m 


ikR 


k:it  R 


X    2    2    2 
i  k  sin  0   sin  (t)  +  — -,  (k  -  k  cos  9 ) 

'      [il  Li. 


U)  1^ 
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(3  6)         f  =  -  ^^    r  e"^^  ^^"  0   sin  (()  _^  i   sin  ([)  -  n  sin  0  ^^-^  ^^^  q  g^^  ^ 
^  Jt  R    L  i   sin  (|)  +  n  sin  9 


kX         X 
where  n  =  — -;  =    j==     may  be  interpreted  as  normalized  surface  reactance, 

Now  we  turn  our  attention  to  the  surface  wave  far  field.    Substituting 
the  last  term  of   (2,l8)   into    (2,10)  and  using   (2,23)^   we  obtain  the   surface 
wave  field; 

"        X  e  "i^    -\Y  +  i    LA^   +  ^^    |x|    +  s   zj 


(3,7)        f.     =^      j  e-^y-^^LV^     +^      |x|+szj      ^^ 

-  k 


'  2^     ^         \^^X2 


X      p  '^        o        ^ 

It  should  be  remembered  that  X  =  —  K     and  K     =  k     -   s  .      In  order 

CO  [I 

to  evaluate  the  integral  by  the  method  of  stationary  phase  it  is  convenient 
to  make  the  following  substitutions 


(3.8)  z     =  p  cos  ii,    |x|    =    P  sin  ii ,     -     =     o. 

(3.7)  becomes 

1 

(3.9)  fg     =     /     A(^,  n,   a)  e^  ^  P     ^^^^  "'   « 


where 


2, 


(3.10)  A(-/.,  n,   a)     .    -^      ■        ''"'  ^^^  "  ^^^    ^^  ^  ^^ 

2it     x/i  -  a2  +  n2(l^a2)2 

and 

(3.11)  P  {^,n,a)  =  a  cos  i;^/     +    1  -  a     +  n^  (l  -  (I'f'         sin  -^ 
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when  k  p  is  sufficiently  lairge 


(3.12)     f^  =     A(a  ) 
s  •= 


2  Jt 


kp 


-i2     . 


i   [k  p  P(^,n,a J"  l] 


It  can  be  easily  shovn  that  ■„  P('^,n,ag)  is  always  negative. 

stationary  point  a       must   satisfy  the  following  equation: 


The 


(3.13) 


dP 
da 


cos  ij/  -   sin  tjj 


Qig  [l  +  2n   (l-Qig)  J 


=  0 


We  may  easily  express  ij/  in  terms  of  n  and  o; 


(3.1^)         tan  Tp 


2  2 

1  -  a^     + 


2   /-,       2  v2] 
n     (I-q;^  )  J 


a 


[l  +  2  n^  (1-a     )  ] 


The  right  hand  side  of  equation  (3.1^)  is  a  monotonically  decreasing 

function  in  the  interval  0  <  a  <  1     for  all  n.     The  truth  of  this  statement 

is  siiggested  hy  Fig.    2  and  can  be  proved  by  differentiating   (3.1^)  with  respect 

to  q;  .      Therefore  a     can  be  found  numerically  for  the  directions  of  observation 
s  s 

0  <  ip  <  90  ,     Those  values  coiresponding  to  the  interval  90     <  '>P  <  I80     and 

-1  <  a     <  0  are  obtained  by  the  symmetry  requirement.      Now  the  surface  wave  far 

field  may  be  written  as 


(3.15)  f    =  -^    ^ 

s        x/2jtkp 


'  2    , 

a    + 
s 


a. 


[l+2n2(l-a^)]  +ri-a2  +  n^(l-a^)^l  ■      .T^,,  .  « 

"-  s -1      L      c;         ^       c:'' jj     i[kpP(^,n,a  )-]7 


l  +  2n^(l-a^)l  ri-aj4  n^(l-aj)^]i 


P 


g-nk    (1-a^)    (y  +h) 
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where  OC     is  a  fimction  of  ip  as  indicated  by   (3.1^).      The  electric   and  magnetic 
field  components  of  the  far-zone   surface  wave  may  be  obtained  by  substituting 

(3.15)  into   {2.h)  and   (2.5).      However ,    in  order  to  understand  the  power  flow 
of  the   surface  wave  far  fields    it  is  desirable  to  use  the  cylindrical  co- 
ordinates   {P,ii,j).      We  first  write  down  the  components  of  f   : 

/  V  ^  ■^  ^ 

(3.16)  f     z     =     f     cos  ^  p   -   f     sin  ijj     ip 

So  s 

Substituting   (3.I6)   into   (2.3)  yields 

(3.17)  E     =  p    |_-n  k(l-a^)   sin^f    J    +    ?//   |_-  n  k   (l-a  )  cos  ^     f    J 

+  y   I    -   i  k  P  sin  ip  f     -   i  k  cos  '^  f    J 


(3.18)     H  =  ^-1 


jk^  P^     sin  ^  f     +  k^   (^  )^  cos  t  f     -  n\^(l-a^)^  cos  f  f  ] 
I-  ^ip  ^      s  ^hip  ^      s  '        s'  ^      s-^ 


^ 


+  ^     n\^    {l-C^f  sin  li/  f     -  k^P^  sin  ?//  f     -  k^P.— ,  cos  ip  f 
^    ^  s  ^s  ^s  d()[/s-' 

^r  22  22SP  ~I 

+  y   I    -   i  n  k   (l-o;   )P  cos  ip  f     +  i  n  k   (l-a   )  —  sin  ■^p  f      \ 
s  s  s      ^1^  s 

1 
where  all  terms  of  higher  order  than  p       have  been  neglected.      The  average 

Poynting  vector  is 

(3.19)     iRe(ExH*)      =     i  Re    (E^  Hy     -   E^  H^    )   p 

+  2  Re    (E     H       -EH      H'   +  2  Re    (E     H,      -   E,    H      )  y 
^yp  py^^        '^  ^p-^  ip     p    '   -^ 
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Lengthy  but  straightforward  substitutions  into  the  above  equation  show  that 
the  f   and  y  components  vanish  and  give  the  exj^ression  for  the  p  component 

—  (l-a  )  )   a 
av  p      CO  i-i      s     s 


I     .2  r  2 

(3.20)    (P   )   =  i  HSL  (i_a2^   a^  r^  ^  2  ^2  ^_^2^  j  _^  r_^2  ^  ^2 ^  2^21 


where  f  is  given  in  (3.15). 
s 


Integrating  (3. 20)  with  respect  to  y  gives  the  pattern  function  for  the 
total  surface  wave  power. 

^    n(l-aj)^  |a^[l+2n^  (1-a^)]  +[l-a^+n^(l-a^)^] 


^3-21)  J    (P^v  ^p  ^y  =8  ,  ,,p        a2  ^  [1  -,  2n2(l-a2)]  n_«2  ^  ,2  _^2  27 
Q  si-       's'-JLs^s-l 


exp    L-2nk     (l  -  a^)h] 
This  is  a  function  of  ip  through  equation   (3.l4). 


k.      Discussions 

Both  the  phase  and  amplitude  of  the  excited  surface  wave  found  in  the 
last  section  are  dependent  upon  the  direction  of  observation.   The  surface  wave 
amplitude  shows  the  typical  exponential  decay  away  from  the  reactive  surface. 
Notice,  however,  that  the  exponential  decay  factor  also  depends  upon  the  di- 
rection of  observation.   Because  the  direction  of  observation  is  approximately 
constant  in  any  fixed  small  region  far  away  from  the  source,  the  three  di- 
mensional surface  wave  may  be  regarded  locally  there  as  an  oblique  plane  surface 
wave.   It  is  recalled  that  none  of  the  anisotropic  characteristics  of  the  dipole- 
excited  surface  wave  are  present  in  the  two  dimensional  treatments  where  an 
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infinite  line  source  is  assumed  '   or  where  the  propagation  of  an  ohlique 

plane  surface  wave  is  considered  .   In  the  case  of  a  finite  line  source  we 

may  expect  to  observe  the  two  dimensional  behavior  at  points  close  to  the 

line  source;  however,  the  three  dimensional  behavior  should  appear  when  the 

distance  from  the  source  is  large  compared  with  the  length  of  the  source. 

There  must  be  a  transition  between  the  two  dimensional  and  three  dimensional 

character! stic  s . 

The  surface  wave  phase  pattern  and  the  surface  wave  power  pattern 
have  been  computed  numerically  for  different  values  of  normalized  surface 
reactance  n  and  h  =  0.   The  phase  curves  using  rectangular  coordinates  are 
plotted  in  Figure  3  for  0  <<;!/<  9O  .   Those  corresponding  to  9O  <  -^  <  I80 
may  be  obtained  by  the  symmetry  requirement.   The  phase  curve  for  n  =  0.1 
is  very  close  to  unity.   The  surface  wave  power  patterns  using  polar  coordinates 
are  plotted  in  Figures  k-   and  5  foi"  a  magnetic  point  dipole  of  unit  stren^h. 
The  surface  wave  power  density  vanishes  at  ^  =  0  and  I8O  for  all  values  of  n. 
It  is  interesting  to  notice  that  the  surface  wave  power  patterns  for  n  =  0.1  and 
0.5  are  broadly  directional  with  the  maximum  at  ^  =  9O  ^  while  those  for  n  =  5 
and  10  are  sharply  directional  with  the  maximiun  at  an  angle  close  to  tjj  =  0 
and  180  .   This  observation  may  lead  to  the  conclusion  that  the  surface  wave 
power  pattern  of  a  finite  line  source  resembles  that  of  a  continuous  array  of 
isotropic  point  sources  for  small  values  of  n.   However,  a  similar  statement 
is  not.  valid  in  the  case  of  large  values  of  n. 

If  frequency  is  very  low,  the  limiting  form  of  the  wave  equation  and 
boundary  condition  in  our  problem  will  be  reduced  to  the  same  as  those  in 
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Kelvin's  ship  wave  problem.  '   But  we  are  interested  in  the  fields  above 
the  reactive  plane,  while  the  ship  wave  solution  deals  with  the  space  below 
•the  water  surface.   Although  the  wave  amplitude  and  phase  in  both  problems 
depend  upon  the  angle  of  observation,  we  find  one  wave  in  all  directions 
but  vanishing  at  ij/  =  0     and  l80  ,  while  the  ship  wave  consists  of  two 
orthogonal  waves  and  is  confined  within  a  sector  behind  the  ship. 
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